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Abstract 

We consider operators of the form H + V where H is the one-dimensional 
harmonic oscillator and V is a zero-order pseudo-differential operator which is 
quasi-periodic in an appropriate sense (one can take V to be multiplication by 
a periodic function for example). It is shown that the eigenvalues of H+V have 
asymptotics of the form \ n (H + V) = X n (H) + W(yfn)rT 1 ^ + 0(rT l l 2 ln(ra)) 
as n — * +00, where W is a quasi-periodic function which can be defined 
explicitly in terms of V. 



1 Introduction 

The one-dimensional harmonic oscillator is the operator 

H = "J? + (axY - 

where a is a positive parameter. We can consider H as an unbounded self-adjoint 
operator acting on L 2 (R). The determination of the spectrum of if is a classical 
problem — virtually any introductory book on quantum mechanics has a section 
devoted to this topic. In particular H has a compact resolvent and hence a discrete 
spectrum. Furthermore, the eigenvalues of H are simple and can be enumerated as 

\ n (H) — a(2n+ 1), n e N . 

A normalised eigenfunction corresponding to X n (H) can be chosen as 

a 1/4 



-ax 2 j2 



where H, n is the n-th Hermite polynomial. 



H n (\/ax) 
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The purpose of this paper is to study the large n asymptotics of the eigenvalues 
of the perturbed operator H + V when V is a self-adjoint quasi-periodic pseudo- 
differential operator of order 0. More precisely, we assume V can be written in the 
form 

aeA 

where A C T*IR = M 2 is a countable discrete index set and, for each a = (a x , a^) G 
T*R, we define U a to be the unitary operator on L 2 (IR) given by 

U a (f)(x) = e ia * a t /2 e ia * x <j)(x + a ( ). (3) 

The Va's are just complex coefficients. 

Since U* = C/_ a for any a G T*R, the condition that V is self-adjoint can be 
rewritten as the requirement 

a G A => -a G A and V- a = 14, a G A. 

We will also assume the V a s satisfy the following condition (essentially a regularity 
assumption) ; 

^|a| 3 |\/ a | < +00. (4) 

aeA 

In particular, this condition ensures that the right hand side of (0) is absolutely 
convergent in operator norm, making V a well defined bounded operator. Since H 
has a compact resolvent the same must then be true for H + V; it follows that the 
spectrum of H + V also consists of discrete eigenvalues. 

Remark. If we take A = {(um, 0) |m £ Z} then V is the operator of multiplication 
by a function with period to whose m-th Fourier coefficient is simply k/^Vj^m.o)- 
Condition (@J) becomes a standard regularity requirement (that the function V should 
be a "bit more" than C 3 ). 

In general we may consider V to be a zero-order pseudo-differential operator 
with Weyl-symbol XlaeA V a e^ axX+a ^ (n.b., U a is the operator with Weyl-symbol 
e i(a x x+a^)y if A is a rational periodic lattice then V will be a periodic operator (in 
the sense that it commutes with a specific translation operator). Taking A to be 
an irrational periodic lattice, or an irregular discrete set, leads to a generalisation 
of such periodic operators; when we apply "quasi-periodic" to V we mean this 
particular type of generalisation. 

If G A then the corresponding term in V is Vq times the identity operator and 
will thus cause a simple shift in the spectrum of H by Vo- This term is included 
in the statement of the main result (Theorem 11.11 below) but thereafter we shall 
assume V = 0. We also set A' = A\{0}; since A is discrete, T*M\A' contains a 
neighbourhood of 0. 
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Define a metric \-\ a on T*R by |a| a = (a 1 a 2 + aa^) 1 ^ 2 . This metric is equivalent 
to the usual metric | • | so condition can be rewritten as 

NalKl < +00 for all p < 3. (5) 

aGA' 

The main result of the paper is the following. 

Theorem 1.1. Suppose V given by (J2J) satisfies (JIJ) (or equivalently ©j. Then the 
eigenvalues of the operator H + V satisfy 

X n (H + V) = a(2n+l) + V + W(^)n~ 1/4 + 0(n- 1/2 ln(n)) 

as n — > 00, where W : 1R — > R is the quasi-periodic function defined by 



W(X) = 




The presence of the quasi-periodic function W means the first order asymptotics 
given by Theorem 11.11 contain considerably more information about the operator V 
than one might expect (c.f. the simple power type asymptotics for the case when V is 
given as multiplication by an element of ([EE]) or for the operator —d 2 /d8 2 +V(6) 
on S* 1 (see Theorem 4.2 in |MU| D. In particular we note that if V is given as 
multiplication by a periodic function, knowledge of the first order asymptotics of 
X n (H + V) allows the Fourier coefficients of V to be "half" determined (the values 
of V(_ ma;i o) + V(mu,o), m E N, can be determined from W). 

It is likely that there exists a full asymptotic expansion for X n (H + V), involving 
further terms with quasi-periodic functions multiplying increasingly negative powers 
of n. Judging by numerical evidence (for example with the potential V(x) = cos(x)) 
the second term in the asymptotics is 0(n _3//4 ). This order (even as an improvement 
of the remainder estimate in Theorem II. lj) appears to involve reasonable subtle 
cancellation effects within the series giving the second term of the asymptotics; no 
attempt to deal with this analysis is made here. 

Remark. With an obvious modification to the definition of W and a remainder 
estimate of 0(n~ 1//3 ln(ra)), Theorem 11.11 also holds for operators V of the form 

V= V a U a d 2 a where K satisfies / (\a\^ 3/2 + \a\l)\V a \ d 2 a < +00. 
Jt*r Jt*r 

In this case V is a pseudo-differential operator of order zero whose Weyl-symbol has 
Fourier transform 2irV a . The | a| ^, term in the condition on V a is then a regularity 
condition, while the |a| a 3 ^ 2 term is a generalisation of quasi-periodicity. 
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The proof of Theorem 11.11 is given in Section 0] using standard ideas to express 
the eigenvalues of H + V in terms of a series involving the resolvent of H and the 
operator V . The non-triviality of Theorem II. II is contained in technical results used 
to establish the convergence of these series. These results are obtained in Sections [21 
and El estimates for the elements (Vtfik, <j>k f ) of the matrix of V with respect to the 
eigenbasis | k G No} are obtained in the former and are then combined to give 
resolvent estimates in the latter. 

Notation. We use C to denote any positive real constant whose exact value is not 
important but which may depend only on the things it is allowed to in a given 
problem. Appropriate function type notation is used in places to make this clearer 
whilst subscripts are added if we need to keep track of the value of a particular 
constant (e.g. Ci(V) etc.). 

We use ||T||, ||T||i and ||T|| 2 to denote the operator, trace class and Hilbert- 
Schmidt norms of the operator T respectively. 



2 Estimates for Matrix Elements 

The aim of this section is to obtain the necessary estimates for the matrix elements 
(V(f>k, 4>k') for all k, k' G N . In turn these will be estimated via 

U k > k ' := {UMM (7) 

defined for all a G T*M and k, k! G No- Since the operator U a is unitary we imme- 
diately get 

\Ut k '\ < 1- (8) 

To obtain more precise estimates we can use the following special function identity 
(see 7.377 on page 844 of [GFO]) to find an explicit formula for U k ' k ; for any < 
k <k! and y, z G C we have 

/ e~ x2 n k {x + y)H kl {x + z)dx = 2 k ' ^k\z k '- k L[ k '- k \-2yz), (9) 
Jm. 

where l^ k ~ k ^ is the generalised Laguerre polynomial. 
Lemma 2.1. For any < k < k! and a G T*M\{0} we have 



U * k ' = \ ^(V2pe*r- k e- p2 L{ k '- k \2p>) 



for some 9 G M, where 

P = 7;{ — + aa i) = ^l a U- 



-l — + aai) = - 



4 



Proof. Introduce the complex number 



y/aac . a x 
uj = — i ■ 



From (J7J), © and (0) we get 

—n-(k+k')/2 r 

ia x a ( /2 j ^ia x x ^-a{x+a^)' 2 /2 ^-ax 2 /2 



7T 



o-(fc+fc')/2 

* — aa 2 /2+ia x a^/2 



l/2+ia x a ( /2 f & -x 2 ^ ^ _ ^ + ) ^ ^ _ ^ ^ 



^ 2 ( fe '- fc V2(_ w )fc'-fc e - 2 --!/2+^^/2 jL (fe'-fe)(_ 2a;(a; _^^ ) ) 

where the last line follows from Q. Now |u;| = p while 

2 2 9 2- 

2 aa^ m^a^ aa^ ia x a^ ta x a^ 2 

^ ~ ~2~ + 2 ~ ~4~~4^~^ 2~ + 2 ~ "'^ 

and 

— 2cj(u; — ^/aa^) = — 2u;(— ZU) = 2|a;| 2 . 
The result follows. ■ 

Throughout the remainder of this section we will assume a e T*M\{0} is fixed 
and p > is given by (fTU|) . 

Laguerre polynomials can be expressed in terms of the confluent hypergeometric 
function; using 22.5.54 in [AS j we get 

L[ k '- k \2p*) = QM(-k,k'-k+l,2 P *). 

The confluent hypergeometric function can, in turn, be written as a pointwise ab- 
solutely convergent series of Bessel functions; from 13.3.7 in jASj we get 

M(-k, k' - k + 1, 2p 2 ) = (k 1 -k)\e p2 (p 2 (k' + k + l))~ (k '~ k)/2 

^{ { k' + k + iy^ J Jv-x+iW + k + i), 

where 

A = l, A 1 = 0, A 2 = ^{k'-k + l) (11) 



and, for j > 2, 

(j + l)A j+1 = (j + k' — k)Aj-i -(k> + k + l)A d - 2 . (12) 
It follows from Lemma f2. II that 

v y = jW-W^F^ ^4 ( {k , + k P + 1)1/2 j J k ^ k+j (2 P V¥+k + T), (13) 



3=0 

where 

k'\/ 2 \k'~ k 
k '' k ''~ ~kl\k' + k + l) 
The next two results give estimates for the constants appearing in ()13j) . 

Lemma 2.2. Suppose k' > 2 and < k! — k < k' 2//3 . T/ien 

|A,-| < {k' + k + l) j/3 . 

Proof. Set m = k' — k and n = fc' + k + 1 so 

< m < A;' 2/3 < (fc' + jfc + I) 2 / 3 = n 2 ' 3 

while k' >2 and > so n > 3. 

We have A = 1 = n°, At = < n 1 / 3 and m, 1 < n 2 / 3 so A 2 = |(m + 1) < n 2 / 3 . 
Now let J > 2 and suppose the result hold for j < J. Since 

. J +m n . 

Aj+i = 4r_, Aj_2 

j+i J+ 1 j 1 J+ 1 j 2 

we then get 

1 J+ 1 - J + l J + l J + l 

Now mrT 2 ! 3 < 1 while 

n > 3 =^ n~ 2/3 < 3~ 2/3 < i 

=> J(l - n' 2/3 ) > 1 (as J > 2) 
=> 1 + Jn" 2/3 < J. 

Thus (J + m)n~ 2 / 3 + 1 < J+ 1. Therefore — n^ J+1 ^ 3 and the result follows 

by induction. ■ 

Lemma 2.3. If < k < k' then F k , tk < 1. 
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Proof. We have 

k'(k' — 1) ... (A; + 1) 



\{k' + k + l)...\(k' + k + l) 



where the numerator and denominator both contain k' — k terms. Now set m 
— k — 1) and n = + A; + 1) so m < n while 

(n + m) (n + m — 1) (n — m — 1) (n — m) 

-Tfc'.A: = • • • • 

n n n n 

Ilk' — k is odd this can be rearranged as 

(n + m)(n — m) (n + m — l)(n — m — 1) n 
n 2 ri 2 n 



Fk 1 



while if k' — k is even we get 

(n + m)(n — m) (n + m — l)(n — m — 1) (n + |)(^ — \) 



k'k 



n 2 n 2 n 2 



The result now follows from the fact that 

(n + m')(n — m!) n 2 — m' 2 



n 2 n 2 



< 1 



for any < m' < n. ■ 

Next we obtain some estimates for the Bessel functions appearing in (fTTTf) . 
Lemma 2.4. For any x, e > and n e [0, x/2] 

Ait e 

{0 e [0,7r] tacos(0) -n| < e} < . 

1 L 1 J 1 3 x 

Proof. Set S = e/x, y = n/x and VL V) s = Cos _1 ([y — S, y + 5]); we need to show that 
\tt y J\ < AnS/3. 

Now set 0o = Gos^ 1 (y) and let £(9) denote the affine function with £(0) = 1 and 
^(0o) — V- It is easy to see that |cos(0) — y\ > \£{0) — y\ which implies \VL Vi s\ < 25/|L| 
where L is the gradient of £{0). On the other hand, y G [0, |] so the minimum value 
for \L\ occurs when y = 1/2; hence 1/\L\ < 2Cos _1 (l/2) = 27r/3 and the result 
follows. ■ 

Lemma 2.5. For any n 6 No and x > 2n we have \J n (x)\ < 4x -1 / 2 . 

Surely this estimate (or an improvement) lies in a book somewhere! 
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Proof. Define a function by f{6) = xsin(#) — n6 so we have the following integral 
representation for the Bessel function J n (see 9.1.21 in jASJ); 



7T 



Jn{x) = ~ / COs(/(0))d0. 



(14) 



Now set 



Q Q = [9 e [0,tt] I \f'(6)\ < x 1/2 } and ^ = [0,7^0 
so J n (x) — (Jo + /i)/7T where = J^ k cos U W) ^ for = 0, 1. Lemma EH gives 

4tt 



M < l^ol < y^ 1/2 - 



(15) 



On the other hand 
h 



sin(/(g)) 
/'(*) 



+ 



sin(/(0))d0. 



0l (/'(*)) 



Now f'(0) = -xsin(0) < on [0,vr] while (/' (0)) 2 > on Thus 



sin(/(0)) d0 



< 



ng) 



d0 



lf'(0) 



an! 



Furthermore f'(9) is decreasing on [0, n] so Qq consists of a single interval. Hence 
<9f2i\{0, 7r} contains at most 2 points. Since /(0) = and /(7r) = — wr we then get 



\h\ < 



sin(/(g)) 
/'(*) 



.9^1 



+ 



/'(*) 



< 6 max . 
an, e ^ \f 



1 



(16) 



Combining (JTSJ), fljft we now get 



1^)1 < -(|/o| + |/i|) < -(^ + 6)x- 1 / 2 < 4^, 

completing the result. ■ 

Lemma 2.6. Suppose k' > 2, < k' - k < p(k' + k + 1) 1/2 and 2p < (k' + k + l) 1/6 . 
Then 

\Ut k '\ < (4(2p)- 1 / 2 + |(2p) 2 )(A; / + A; + l)- 1 / 4 . 
Before starting, note that as a clear consequence of (|14|) we have 



|J n (x)| < 1. 



(17) 
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Proof. Since 2, k < k' 

k'-k < ^k'+k+i) 2 / 3 < \(iy /3 k^ 3 < k a /\ 

Now combining (J13|) with (jllj) . (J 17)) and Lemmas 12.21 and 12.31 we get 

OO j 

\ut k '\ < v^kl E 1^1 (y + T+ iy/2 \Jk'- k+j (2 P Vk' Tk + T) | 

OO 

< |j fc ,_ fc (2pVfc' + fc + 1) | + (fc' + fc + l)" i/6 

< \j k ^ k (2 P VFTkTT)\ + l(2p) 2 (k' + k + iy 1/3 , 

where the last line follows from the hypothesis that p(k' + k + l) -1 / 6 < 1/2. Lemma 
12.51 can now be used to estimate the remaining Bessel function term. ■ 



Main estimate 

The next result is the main estimate we will need for the matrix elements | (V (/))., <pk')\- 
This estimate is valid in a parabolic region around the diagonal k = k'; the width 
of this region is governed by the quantity 

7 := min laL , 

aeA' 

which is positive since A' is discrete and doesn't contain 0. Although not required 
in this paper, we remark that for a general parabolic region around the diagonal one 
is restricted to estimates of the form \ (V(j)k, <f>k')\ < C(V)(k' + k + l) -1 / 6 . 

Proposition 2.7. Suppose V satisfies condition (JHJ) and set 

k = min{l/3,7/(2v / 3)}. (18) 

If n EN and k, k' G No satisfy \k — n\,\k' — n\ < nn 1 ^ 2 then 

\(V<f> k ,<f> v )\ < C{V)n'^. (19) 

Proof. We have \{V(fik, (fik')\ < ||^|| for any k,k' G No so we can increase C(V) if 
necessary to ensure that ()19|) is satisfied for n = 1,2. Furthermore V is self-adjoint 
so | (V4>k>, 4>k) \ — I (V(pk, 4>k')\- It thus suffices to prove the result assuming n > 3 and 
k', k G No satisfy k' > k and \k — n\, \k' — n\ < nn 1 ^ 2 . Then k',k > n — |n 1//2 > |n 
so k' > 2, 

4 

k + k + l > -n (20) 
9 



and 

< k'-k < 2/m 1 / 2 < l(k' + k + 1) 1/2 . (21) 
Now set K ={k' + k+ l) 1 / 6 . Using ©, © and © we have 

\{V4>M\ < Ei^'ilKI < £ 1^11^1+ £ IKI- (22) 

aeA' aeA' aeA' 

|a| a <AT |a| Q >A- 



Since 1 < i\:- 3/2 |a|a /2 whenever | a| a > K, © and (j20J) give us 

E iki < k-^ i a i' /2 i^i < cw^ 1 



/4 



aeA' aeA' 
|aL>if 



Now let a 6 A'. Since |a| Q = 2p (see ()1U|)) the definition of 7 implies 7/2 < p and 
thus k! -k< pK 3 by fl2U). Lemma © and then give 

aeA' aeA' 
|a| a <.K" 

The result follows. ■ 
First order term 

The next result is used to obtain the explicit form for the first order correction term 
in the asymptotics for X n (H + V). 

Proposition 2.8. Suppose V satisfies condition ©. Then 

(V<f> n ,<f>n) = WiMn-^ + Oin- 1 ' 2 ) 
as n — > +00, where W is defined by ©. 

Proof. Let a G T*1R\{0} and set p = |a| a /2. Using (JEU) and the fact that F„ >n = 1 
we get 

u » n = £ % n + 1V / 2 ^(W^). 

Now suppose 2p < N where N := (2n + l) 1 / 6 . Using (fTT|) . (JTTj) and Lemma l2~2l we 
have 

00 

\U2' n - J (2pv / 2^Tl)| < |p 2 (2n + l)- 1 + E/> J '(2^ + l)^ /6 

i>3 

< ±|a| 2 (2n + l)" 1 + i|a|3(2n + l)- 1 /2. 
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Standard asymptotic forms for Bessel functions (see 9.2.1 in |ASj ) give us 

Y 



Mz) 




cos (z J + 0{z~ 



7TZ V 4 



■3/2^ 



while 



d ( 1 / 7T 

— — ■= cos 2; - — 
dz\^z V 4 



2 



and 2pV2n + 1 - 2p v / 2^ < 2~ 1 / 2 pn- 1 / 2 . It follows that 

J (2pV2n + 1) - W-(2p)" 1/2 (2n)" 1/4 cosf2pv / 2^- ^ 
V vr V 4 



< C(2p)- 3 / 2 (2n + I)- 3 / 4 



+ 



'(2p)- 1 / 2 (2n)- 1 / 4 + ^(2pr 3 / 2 (2n)- 3 / 4 ) 2~ 1 / 2 pn- 1 / 2 
< C((2p)- 3 / 2 + (2p) 1 / 2 )n- 3 / 4 . 

Combining the above estimates we thus obtain 



jjn,n 



2 l/4 



-1/2--1/4 



7T 



n 7 cos a Lv2n 



< ^(|a|- 3/2 +|a|> 



whenever |a| a < AT. Using (J2J), ©, © and (jHJ) we thus have 

|(V0 n ,0„)-w(V^)n- 1/4 | 

< Cn-V* (l^ 3/2 + l^)l^|+ £ (l + |a£ 1/a )|K|. 



aeA' 
|a| Q <V 



aeA' 

|aL>7V 



Since 1 < A" 3 |a| 3 < n 1 ^ 2 |a| 3 whenever |a| a > A" the term inside the last sum can 
be replaced with n~ 1//2 (|a| 3 + |a|a )|T4|. Using (J3J) we then get 



\{V4>nAn)-W{Vn)n-^\ < Cn- 1 / 2 E (|a|- 3 / 2 +|a|»)|V a | < C(V)n~ x l\ 
completing the result. 



aeA' 



3 Resolvent Estimates 

For any A G C\a(H) let R(X) — (H — A) -1 denote the resolvent of the operator H; 
we will also write R for R(X) where this should not cause confusion. 
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Let k denote the constant denned in (|18|) . For a given n G N we will make 
repeated use of the partition of N defined by 

I = [k G N | \k-n\ < nn 1/2 } and J = N \I. (23) 

For any e G (0, a) and n G No, let r £in be the anti-clockwise circular contour in 

C centred at A n = X n (H) = a{2n + 1). If A G r £i „ then A = a(2n + l) + ee ie 

for some 9 G [0, 2ir). It follows that |A - X k \ = \2a(n-k) + ee ie \ for any k G N . 
Straightforward arguments then lead to the following estimates; 

^lA-Afcl" 1 < C(e) ln(n), (24) 
k&i 

\ X ~ X k\~ 2 < C( £ )> ( 25 ) 

5^|A- A fc |- 2 < Cn- 1 ' 2 (26) 
fceJ 

and 

|A-A fc |>Cn 1/2 for any k G J. (27) 



The first two results in this section relate to the operator R(X)VR(X), which is 
clearly bounded whenever A is in the resolvent set of H. We show that it is in fact 
trace class while its operator norm decreases as n -1 / 4 for A G T £ ^ n . 

Lemma 3.1. For any n G N and A G r £) „ we have 

\\R(X)VR(X)\\ < \\R(X)VR(X)\\ 2 < C^e)^ 1 ^. 
We remark that since {<ftk \ k G No} is an orthonormal basis of L 2 (R) 

J2 KVfaAk*}] 2 = WVM 2 < \\v\\ 2 . (28) 

fc'eNo 

Proof. Using the orthonormal basis {4>k \ k G No} we have 

\\rvr\\1 = £ \(RVR4>k,M\ 3 = £ i a l -tm?-\\* - m 

A:, fc'eNo fc,fc'GN ' h I I fc ' I 

We will split this sum using the partition (|2*3*|) . Firstly Proposition 12.71 and (|23|) 
imply 

£ c( ">- /2 (Sf^) 2 * 

12 



Now using fl27j), (BED and J2HJ we get 
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^ |A fc - A| 2 |A fc / - A| 2 ^ C " Lp^AjS^IW*'^ 

fceN 1 11 1 fceN 1 1 fc'eJ 
fc'eJ 

< C(y,e)7C x . 



The remaining part of the sum on the right hand side of involves k G J and 
k' G / C No; thus we can estimate this part using an argument similar to the last 
one with k and k' swapped. ■ 

Lemma 3.2. For any n G No and A G r £>n t/ie operator R(X)VR(X) is trace class. 
Furthermore \\R(X)V R(X)\\i is uniformly bounded (in n and X G r e>7 J. 

Proof. The set {0^ I & £ No} is an orthonormal eigenbasis for R with corresponding 
eigenvalues (A& — A) -1 , k G N so implies 

\\R\\l = ^|A fc -A|" 2 < C(e). 

fceNo 

Thus HW^Hi = ||Vi2 2 ||i < ||V|| ||i? 2 ||i < ll-Rlll < C(e)\\V\\. U 

Suppose n G N and j G N. From the previous result we know that R(X)VR(X) 
is trace class for any A G r en . On the other hand R(X)V is bounded (in fact 
\\R(X)V\\ < It follows that 

(R(X)V) j R(X) = (i?(A)V)^ 1 J R(A)V r /?(A) 

is also trace class with trace norm uniformly bounded for A G T £>n . The work in the 
remainder of this section leads to Proposition 13.51 where we obtain an estimate for 
the trace of an integral of such operators. 

Lemma 3.3. Let n > 2, X G r £i „ and suppose f : No — > C satisfies 

\f( k )\ 2 ^ C i and ^ C x n- 1 ^ when kel (30) 

for some constant C\. For each k G No set 

f(k') (Vfafa) 
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(*> = E 



A — Xui 

fc'GNo 



Then there exists a constant K = K(V, e) such that 

l^( fc )| 2 - CfK 2 n~ 1/2 ln 2 (n) and \g(k)\ < dKn- 1 ' 2 ln(n) w/ien kel. 

fceN„ 
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Proof. Since 



E l«K*)l s 

feeN 



E 

fc,fc',fc"eN 



f(k')f(k")(V<f> k r,<f> k )(<f> k ,V(f> k ») 



and 



E (V(p k ',(pk)((pk, V<p k n) 
fee No 



(A — Ajfc')(A — Afc") 
1(^,^)1 < 



it follows that 



E i^)i 2 ^ 



v |/(fc)| v |/(*)| ^ 



kei 



|A — A fc | 



Using the second part of (|3Uj) and (J23j) we get 



E 

fee/ 



l/WI 
|A- A fc 



< Cin-^El^-Afcl- 1 < Ci^^n-^lnCn). 



fee/ 



On the other hand the first part of (}30j) and ()26|> give 



E 

fceJ 



l/WI 
|A-Afe 



< 



(ei/wi 2 ) 1/2 (ei^ 



A; 



1-2 



v fceJ 7 x fceJ 

< CiC 3 n" 1/4 < 2C 1 C 3 n" 1/4 ln(n). 



1/2 



Putting these estimates together now leads to 

£>(A0| a < ClKln^ 
fceNo 



with l£i = || V||(C 2 (e) + 2C 3 ). Now suppose k E I and write = gi(k) + gj(k) 
where 



= E 



/(fcQ (^>fe) 
A — At' 



and flfj(A;) = E 



A — At/ 



fe'e-f fc'eJ 
From Proposition 12. 7\ the second part of (}30|) and (|24jl we get 

|<7/(fc)| < CiC^n-^ElA-^r 1 < C^G^rT^ln^). 



fe'e/ 



1/2 



On the other hand (|2*7j). the first part of (|3Uj) and (j2Sj) give us 

i^(*)i < c(s)n-^(j2\f(n 2 Y\j2\^^M 2 ) 

^k'eJ ' ^fc'eJ ' 

< C 1 G^eW\\n 11 ' 1 < 2C 1 C 5 (e)\\V\\n- 1 / 2 Hn). 

Putting these estimates together now leads to \g(k)\ < CxK^rT^ 2 ln(n) with K 2 
C 4 (V,e) + 2C 5 (e)||V||. Taking K = max{K 1} K 2 } } completes the result. 



14 



Taking f{k) = (V4> n ,(p k ) we can use (J28)l and Proposition 12. 71 to check that (|3(jp 
is satisfied. The next result then follows from Lemma 13.31 by use of induction; we 
can take K = max{||V||, C(V), K'} where C(V) and K' are the constants coming 
from Proposition 12.71 and Lemma f3.3l respectively. 

Lemma 3.4. Suppose n > 2 and j G No- Then there exists a constant K = K(V,e) 
such that for all A G T £ n we have 



E 



(V(f) n ,(f) kl )(V4> kl ,4> k2 ) . . . (V(p k ,4> n ) 



(X-X kl )...(X-X kj ) 



ki,...,kj£ 

Proposition 3.5. Suppose n > 2 and j G N. Then 
1 



n). 



Tr 



2m 



XR(X)(VR(X)y dX 



n 



where K is the constant from LemmaFTJ 



Proof. Since {(p k > \ k' G No} is an orthonormal basis of L 2 (M) we have 

fc'eNo fc'GNo 

Continuing by induction we get 



Xi. — A 



E 



k 1 ,...,k i eNo 



(V(j) k ,(f) kl )(V(j) kl ,(f) k2 ) . . . {Vfa^faj) 
(A fc -A)(A fel -A)...(A fc -A) 



Together with the fact that (R(fi kj , (fik ) = S kjtko (X ko — A) 1 we now get 

Tr(R(VRy) = {R{VR) j K, K) 

fe 6N 

\- (V r feo ,0 fcl )(V A fcl ,0fc 2 ) . . . (V(j) k (j) kj ) 

2^ T\ ~\ n ^ Wfcj.^fco) 



fco,---,fcj6No 



E 



(Afc — A)(A / i Cl — A) . . . (Afcj^j — A) 

(v ^kpAki) (v (j>k 2 ) ■ ■ ■ (v^-iAko) 

(X k0 - A) 2 (A fcl - A) . . . (A^ - A) 



£ A^' 



1=0 



fc ,...,fc J _ieN 
where ^4(A) is the meromorphic function 

(V(p ko ,(i) kl ){V(i) kl ,(j) k2 ) . . . fco ) 



J k0 ,...^m, (^o-A)(A fel -A)...(A fej „ 1 -A) 



(31) 
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Since 

d / J,_1 1 \ J 1 1 17-1 1 1 

dA A (n x^x) = II A^a(II ^Ta 

for any k , . . . , G N , we can rewrite the above equation as 

Tr XR{VR) j = ^(AA(A)) -A(X). 

ClX 

Integrating around the contour r e n it follows that 

1 A{\) dX. (32) 



Tr 2^f XRiVR)3dX = ~hf 



,ri 

The poles of the meromorphic function A(X) occur at the points A = A^ for k G No- 
Since the only such point enclosed by the contour T en is A = A ra , it follows that 
the only terms in the series (|31j) which contribute to the right hand side of (|32j) are 
those with at least one of ko, . . . , kj-x equal to n. With the help of symmetry we 
then obtain the identity 



Tr 2tt? T XR ( VR Y dX 



L a,. _u a, u dX - ( 33 ) 



2m A ™ ~ A fcl ,...^ ie N„ ^ " A ) ' ' ' ^ ~ A ) 



For any A G r £j „ we have |A n — A| — e while 

(V4> n ,(f) kl )(V(f) kl ,(j) k2 ) . . . (V(f> k . i; n ) 



E 



/, i /, (A fcl - A) . . . (A fe ._ x - A) 



< K j n- j/4 ln j ' 



n 



by Lemma f3. 41 Since the length of r e n is 27re we finally get 
Tr — / XR(VR) j dX <— I - K j n~^ 4 ln J ' -1 (n) dX = K j n~ j ^ 4 ln J ' _1 (n), 



completing the result. 

Taking j = 1 in ()33)1 leads to the formula 

Tr / Ai?(A)Ki?(A)dA 



2h 



1 ; 1 (V0 n ,0 n )dA = (V<f> n ,<f> n ). (34) 



27T2 An A„ — A 

This is needed to obtain the first order correction term in Theorem ll.il 
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4 Proof of Theorem 11.11 

Lemmas O and H21 give us \\R(X)VR(X)\\ < Cm' 1 / 4 and ||i2(A)Vi2(A)||i < C 2 for 
all n G N and A G T e>n . In particular || (VR(X)) 2 \\ < \\ VWdn' 1 ^. We also note that 
||i?(A)|| = e^ 1 . It follows that for any j G No we get 

||(ra(A))*|| < ||(W?(A)) 2 |K < (\\VWdn-VV 

and 

||(ra(A))* +1 || < ||v||||i?(A)||||(\/i?(A))^|| < iiyiie-^nyiicm- 1 /^. 

Choose iV'eNso that ||K||C 1 A^'~ 1/4 < 1/2. It follows that for any n > N' and 
A G T e n the series 

oo 

(J + V^A))- 1 = ^(-yitl(A)y (35) 

3=0 

is absolutely convergent and has norm bounded by 2(1 + ||V||£ -1 ). In particular, 
I+VR(X) is invertible with a uniformly bounded inverse for all n > N' and A G T £ n . 
On the other hand, the series 

oo oo 

T(A) := "£R(\)(-VR(\)y = -R(\)VR(X)J2(-VRW) j 

3=1 3=0 

is convergent in trace class with 

ll^(A)iu < n^y^HiiKJ + v^A))- 1 !! < 2c 2 (i + ||y||,- 1 ) 

for n > N' and A G r £jn . Setting 

T n = -^-j AT(A) d\ 
it follows that we have an absolutely convergent expansion 

oo 

TrT n = -VTr — : d> \R(X)(-VR{X)) j dX 

. ^ 2"7f 2 / p 

whenever n > N'. 

Now choose N > N' so that KN -1 ^ ln(iV) < 1/2 where K is the constant given 
by Proposition 13.51 Using this Proposition and the above results it follows that 



oo 

V Tr d> XR(X)(-VR(X)y dX 

j=2 J 1 e,n 



< 2K 2 n- 1 ' 2 ln(n) 
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for all n > N. Therefore 

TrT n = Tr — / XR(X)VR(X) dX + 0(n~ 1/2 ln(n)) 

= (V<p n ,<Pn) + 0{n^ 2 hx{n)) 

for all n > N, where we have used (p?4"j) . 

The argument can be tied together using a standard resolvent expansion. Set 
R V (X) = (H + V - X)- 1 and let n > N. Then 

oo 

R V (X) = J R(A)(l + Ki?(A))- 1 = R(X) (-VR(X)y. 

3=0 

The right hand side of (pI5j) will still converge if V is replaced with gV for some 
g G [0, 1]. Hence a(H + gV) n T £ n = 0. Since the eigenvalues of H + depend 
continuously on g, it follows that T e n must enclose X n (H + 1/) but no other points 
of a(H + K). Thus we can write 

X n (H + V)- X n (H) = Tr I X(R V (X) - R(X)) dX 

„ oo 

= -— Tr * a y2 R ( x )(~ VR ( x )y dx 

u i e,n J = l 

= TrT n = (^0 n ,0 n ) + O(n- 1 / 2 ln(n)). 
Theorem 11.11 now follows from Proposition 12.81 
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